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The B — > X s 77 decay: 
NLL QCD contribution of the Electromagnetic Dipole operator C 7 
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We calculate the set of 0(a s ) corrections to the double differential decay width dr77/(dsi ds2) 
for the process B — > X a j-y originating from diagrams involving the electromagnetic dipole operator 
O7. The kinematical variables Si and S2 are defined as s t = (pt — qt) 2 /m\, where Pb, ?i, 92 are 
the momenta of fe-quark and two photons. While the (renormalized) virtual corrections are worked 
out exactly for a certain range of si and S2, we retain in the gluon bremsstrahlung process only 
the leading power w.r.t. the (normalized) hadronic mass S3 = (pb — qi — q?) 2 /ml in the under- 
lying triple differential decay width dT 77 / ' (dsids2dss) . The double differential decay width, based 
on this approximation, is free of infrared- and collinear singularities when combining virtual- and 
bremsstrahlung corrections. The corresponding results are obtained analytically. When retaining 
all powers in S3, the sum of virtual- and bremstrahlung corrections contains uncanceled 1/e singu- 
larities (which are due to collinear photon emission from the s-quark) and other concepts, which 
go beyond perturbation theory, like parton fragmentation functions of a quark or a gluon into a 
photon, are needed which is beyond the scope of our paper. 



Inclusive rare i?-meson decays are known to be a unique source of indirect information about physics at scales 
of several hundred GeV. In the Standard Model (SM) all these processes proceed through loop diagrams and thus 
are relatively suppressed. In the extensions of the SM the contributions stemming from the diagrams with "new" 
particles in the loops can be comparable or even larger than the contribution from the SM. Thus getting experimental 
information on rare decays puts strong constraints on the extensions of the SM or can even lead to a disagreement 
with the SM predictions, providing evidence for some "new physics" . 

To make a rigorous comparison between experiment and theory, precise SM calculations for the (differential) decay 
rates are mandatory. While the branching ratios for B — > X s ~f [ij and B — > X s £ + £~ are known today even to 
next-to-next-to- leading logarithmic (NNLL) precision (for reviews, see 0,01), other branching ratios, like the one for 
B —> -X" s 77 discussed in these proceedings, has been calculated before to leading logarithmic (LL) precision in the 
SM by several groups 0-01 an d only recently a first step towards next-to-leading-logarithmic (NLL) precision was 
presented by us in [8(. In contrast to B — > X s j, the current-current operator O2 has a non- vanishing matrix element 
for b S77 at order a® precision, leading to an interesting interference pattern with the contributions associated 
with the electromagnetic dipole operator O7 already at LL precision. As a consequence, potential new physics should 
be clearly visible not only in the total branching ratio, but also in the differential distributions. 

As the process B — > X s j~f is expected to be measured at the planned Super B-factories in Japan and Italy, it is 
necessary to calculate the differential distributions to NLL precision in the SM, in order to fully exploit its potential 
concerning new physics. The starting point of our calculation is the effective Hamiltonian, obtained by integrating 
out the heavy particles in the SM, leading to 



I. 



INTRODUCTION 




(1) 
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where we use the operator basis introduced in 

01 = (sL^T a c L ) (c Ll ^T a b L ) , 2 = (s LW ) (c Ll »b L ) , 

O s = (s-L7A)E 9 (97 M 9), 04 = (sL^T a b L )J2 q (qrTaq), 

5 - (sL^ v i P b L )T, q (q^ri p g) , 06 = (sLW»j P T a b L )j: q (<irryT a q) , 

07 = mtifj) [s L a^b R ) , 8 = $p rh b (n) (s L a^T a b R ) G%, . 

The symbols T a (a = 1,8) denote the SU(3) color generators; g s and e, the strong and electromagnetic coupling 
constants. In eq. (j2]), m&(//) is the running 6-quark mass in the MS-scheme at the renormalization scale As we 
are not interested in CP-violation effects in the present paper, we made use of the hierarchy V u bV* s <C VtbV t * when 
writing eq. ([IJ. Furthermore, we also put m s = 0. 

While the Wilson coefficients Ci (fi) appearing in eq. (JTJ are known to sufficient precision at the low scale fi ~ mb 
since a long time (see e.g. the reviews and references therein), the matrix elements (sjj\Oi\b) and (sjj g\Oi\b) , 
which in a NLL calculation are needed to order g 2 and g s , respectively, are not known yet. To calculate the (Oi, Oj)- 
interference contributions to the differential distributions at order a s is in many respects of similar complexity as the 
calculation of the photon energy spectrum in B — > X s ^j at order a 2 needed for the NNLL computation. As a first step 
in this NLL enterprise, we derived in our paper [sj, the 0(a s ) corrections to the (Or, 07)-interference contribution 
to the double differential decay width dT/(dsids2) at the partonic level. The variables s\ and s 2 are defined as 
s i = (Pb — Qi) 2 / m bi where pb and qi denote the four-momenta of the &-quark and the two photons, respectively. 

At order a s there are contributions to dTm / (dsidsz) with three particles (s-quark and two photons) and four 
particles (s-quark, two photons and a gluon) in the final state. These contributions correspond to specific cuts of the 
6-quarks self-energy at order a 2 x a s , involving twice the operator Oj. As there arc additional cuts, which contain for 
example only one photon, our observable cannot be obtained using the optical theorem, i.e., by taking the absorptive 
part of the &-quark self-energy at three-loop. We therefore calculate the mentioned contributions with three and four 
particles in the final state individually. 

We work out the QCD corrections to the double differential decay width in the kinematical range 

0<si<l ; < S2 < 1 - si . 

Concerning the virtual corrections, all singularities (after ultra-violet renormalization) are due to soft gluon 
exchanges and/or collinear gluon exchanges involving the s-quark. Concerning the bremsstrahlung corrections 
(restricted to the same range of si and S2), there arc in addition kinematical situations where collinear photons 
are emitted from the s-quark. The corresponding singularities are not canceled when combined with the virtual 
corrections. We found, however, that there are no singularities associated with collinear photon emission in the 
double differential decay width when only retaining the leading power w.r.t to the (normalized) hadronic mass 
S3 = (pb — qi — q2) 2 /m 2 in the underlying triple differential distribution dT^j / (dsids2ds 3 ). Our results of our paper 
are obtained within this "approximation". When going beyond, other concepts which go beyond perturbation theory, 
like parton fragmentation functions of a quark or a gluon into a photon, are needed [10( . 

II. LEADING ORDER AND FINAL RESULTS FOR THE DECAY WIDTH 

In d = 4 dimensions, the leading-order spectrum (in our restricted phase-space) is given by 
rfTff _ a 2 m 2 ^) m l\C 7 , eS} ^)\ 2 G 2 F \V t bV; s \ 2 Q 2 (l- Sl -s 2 ) 



ds 1 ds 2 1024 7T 5 (1 - si) 2 si(l - s 2 ) 2 s 2 

where 



ro ■ (3) 



r a = -48s^s? + 96s^s?-56s 2 s? + 8s? + 96s^s?-192s^s 2 + 112s 2 s?-56sisi + 
lUsjsi - 96s 2 si + 8si + 8s^ + 8s 2 
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The complete order a s correction to the double differential decay width dYn/(dsi ds2) is obtained by adding the 
renormalized virtual corrections and the bremsstrahlung corrections. Explicitly we obtain 



rfr 



(i) 

77 



a 2 mt(riml\Cr, eff (ri\ 2 G%\V tb V t * s \ 2 Q 2 d ^ a 



dsi ds2 



1024 rr 5 



-4r (l- Sl -.s 2 ) ji_ + ■ 

(1 - si) 2 si (1 - s 2 ) 2 s 2 m b 



(4) 



where / can be found explicitly in 

The order a s correction dT^ /(dsids^) in Eq. ([J} to the double differential decay width for b — > X s "fj was the 
main result of our paper 



III. SOME NUMERICAL ILLUSTRATIONS 



In our procedure the NLL corrections have three sources: (a) a s corrections to the Wilson coefficient CV )e //(ju), (b) 
expressing fhb{p) in terms of the pole mass mj, and (c) virtual- and real- order a s corrections to the matrix elements. 
To illustrate the effect of source (c) , which is worked out for the first time in our paper Q , we show in Fig. Q] (by 
the long-dashed line) the (partial) NLL result in which source (c) is switched off. We conclude that the effect (c) is 
roughly of equal importance as the combined effects of (a) and (b). 

For completeness we show in this figure (by the dotted line) also the result when QCD is completely switched off, 
which amounts to put fj, = mw in the LL result. 

From Fig. Q] we see that the NLL results are substantially smaller (typically by 50% or slightly more) than those 
at LL precision, which is also the case when choosing other values for s 2 . 

In the numerical discussion above, we have systematically converted the running 6-quark mass rat(/x) in terms of 
the pole mass nib- As perturbative expansions often behave better when expressed in terms of the running mass, we 
also studied the results obtained when systematically converting m& in terms of m&(/i). After doing also this version, 
we observe the following: Generally speaking, NLL corrections are not small for both cases, when taking into account 
the full range of /i, i.e., mi,/ 2 < /i < 2 m,b- More precisely, in the MS version they are large for /i = nib/ 2 and smaller 
for larger values of fi, while in the pole mass version they are large for all values of /i. 

We stress that the numerically important contributions involving the operator 2 are not discussed in our paper. 
Therefore, the issue concerning the reduction of the /i dependence at NLL precision cannot be addressed at this level. 
Finally, the relevant input parameters that we used in our analysis together with the values of the Wilson coefficient 
C7 and the strong coupling a s at different values of the scale /i are listed in Table I. 



Parameter 


Value 


m b (polc) 


4.8 GeV 


rat (p o1c ) 


175 GeV 


M W 


80.4 GeV 


M Z 


91.19 GeV 


G F 


1.16637 X 10~ 5 GcV~ 2 




0.04 


a- 1 


137 


c 3 (M z ) 


0.119 







a?,.,., 00 


C 7.ff"" 


M = M w 


0.1213 


-0.1957 


-2.3835 


f = 2 m b 


0.1818 


-0.2796 


-0.1788 


p. = m b 


0.2175 


-0.3142 


0.4728 


l> = m b /2 


0.2714 


-0.3556 


1.0794 



TABLE I: Left: Relevant input parameters . Right: a s (fi) and the Wilson coefficient Cj^ff(ii) at different values of the scale fi. 



IV. CONCLUDING REMARKS 

In the present work we calculated the set of the 0(a s ) corrections to the decay process B — > ^77 originating 
from diagrams involving O7. To perform this calculation, it is necessary to work out diagrams with three particles 
(s-quark and two photons) and four particles (s-quark, two photons and a gluon) in the final state. From the 
technical point of view, the calculation was made possible by the use of the Laporta Algorithm to identify the needed 



FPCP2012-62 



Presented at Flavor Physics and CP Violation (FPCP 2012), Hefei, China, May 21-25, 2012 




FIG. 1: Frames l)-3): Double differential decay width 0^77 / (dsid,S2) as a function of si for S2 fixed at S2 = 0.2. The dotted(black), 
the short-dashed (red) and the solid line(blue) shows the result when neglecting QCD-effects, the LL and the NLL result, respectively. 
The long-dashed line(purple) represents the (partial) NLL result in which the virtual- and bremsstrahlung corrections worked out in our 
paper [8| are switched off. In the frames 1), 2) and 3) the renormalization scale is chosen to be fj, = mi/2, fi = mi and fj, = 2mj, 
respectively. Down right: The relevant phase-space region for s\ and S2. 

master integrals and by applying the differential equation method to solve the master integrals. When calculating 
the bremsstrahlung corrections, we take into account only terms proportional to the leading power of the hadronic 
mass. We find that the infrared and collinear singularities cancel when combining the above mentioned approximated 
version of bremsstrahlung corrections with the virtual corrections. The numerical impact of the NLL corrections is 
large: for dTn/(dsi ds2) the NLL result is approximately 50% smaller than the LL prediction. 
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